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Abstract 



Using photon-photon collisions as a particularly transparent study case we investigate the production 
mechanisms for Be mesons. In nonrelativistic approximation and to 0{a'^a1) it is shown that recom- 
bination of h- and c-quarks dominates by far over h and c fragmentation. This dominance persists 
up to the highest accessible transverse momenta and leads to distributions in energy which differ 
completely from the spectra expected on the basis of the fragmentation functions. For processes in 
which a 66-pair is radiated from a primary c-quark, the fragmentation description is found to be 
inadequate. We anticipate important implications of these results for hadronic production of heavy 
quark resonances. Using realistic photon spectra we predict two-photon production rates for Be and 
B* at present and future e"'"e~ -machines. 
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1 Introduction 



Since the top quark seems to be too short-lived for quarkonium-hke resonances to form, the 
6c-bound states are most hkely the only narrow heavy quark resonances besides the well-known 
charmonium and bottomonium systems. Till now these so-called mesons have not yet been 
observed. However, they are expected to be in experimental reach with some luck already at 
LEP, more likely at the TEVATRON, but most certainly at the LHC In addition, also the 
HERA-B experiment may get a glimpse on these very interesting bound states. 

Because of flavour conservation in strong and electromagnetic processes, the Be ground state 
can only decay weakly unlike the r]c and rjb. This fact provides unique possibilities to investigate 
important aspects of strong and weak interactions and their interplay. Indeed, the nonrela- 
tivistic nature of these systems allows to calculate also genuinely nonperturbative quantities 
such as total production cross sections, fragmentation functions, decay constants, and transi- 
tion form factors. For this reason, the Be system can serve as a testing ground for quantitative 
approaches to confinement problems and our understanding of soft hadron physics. 

So far, the estimates only include the contributions of lowest order in the relative momentum 
of the constituent quarks and in the strong coupling constant. Moreover, they are subject to 
numerical uncertainties in the bound state wave functions, relevant quark masses and coupling 
constants. Independently of that, the calculations have not yet reached complete agreement 
on the hadronic production cross sections. Also, there are different opinions on the relative 
importance of recombination and fragmentation mechanisms, such as be —>■ Be and b —>■ BcC, 
respectively. 

In this paper, we take up the latter question and clarify the issue of recombination versus 
fragmentation. In order to liberate our study from unnecessary complications, we consider 
photon-photon instead of gluon-gluon scattering which is the dominant subprocess in hadronic 
collisions. The two reactions have many features in common including the competition of 
recombination and fragmentation processes. On the other hand, the absence of coloured 
quanta in the initial state of 77-processes drastically reduces the number of Feynman diagrams 
at a given order and makes the calculations much more transparent. We realize, however, that 
the nontrivial colour coefficients in (7 (yf-pro cesses somewhat change the relative weight of the 
different mechanisms in comparison to 77-scattering. 

Preliminary results of our study were already reported in ref. However, there the focus was 
more on the prospects of observing Be mesons in e"'"e~-annihilation. Here, we complete our 
analysis of Be production in 77 collisions and also correct a normalization error in the fragmen- 
tation estimate given in ref. Furthermore, we show the relevant differential distributions 
and point out the features which distinguish recombination from fragmentation. In addition, 
we examine in detail the theoretical validity of the factorized description in terms of heavy 
quark production, 77 —>■ cc and bb, and subsequent fragmentation, c — > Beb and b —>■ BcC, 
respectively. Finally, we present predictions on the total cross section for Be and B* produc- 
tion in collisions of high energy photon beams which may be obtained by back-scattering of 
low energy, high intensity laser beams at future Linear Colliders, and in collisions of softer 
bremsstrahlung photons. 
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When our work was near to completion, a related study appeared in ref. 0. If we use the 
same parameters we reproduce the total cross sections within the Monte Carlo uncertainties, 
except for the highest 77 cm. energy of 100 GeV considered in ref. where we disagree by 
ten per cent. Also the differential distributions presented in ref. [0 show in general the same 
behaviour as ours. 



2 Classification of processes and calculation 

In 77 collisions. Be mesons are dominantly produced in association with a b- and c-quark jet, 

7(^1) + 7(^2) B,{P) + b{p,) + c(pj. 

Here, the particle momenta are indicated in parentheses. To lowest order in the electromag- 
netic and strong coupling constants a and a^, respectively, the above process is described by 
twenty Feynman diagrams. These can be classified in different topologies as indicated in Fig. 1. 
The diagrams belonging to set (I) can further be distinguished by the flavour of the primary 
quark line coupled to the photons. This line carries either the bottom (1^) or the charm (I^) 
flavour. The complete subsets of six diagrams with primary b- and c-quarks are found by 
interchanging the initial photons as required by Bose symmetry. These subsets are individu- 
ally gauge invariant. The remaining eight diagrams belonging to set (II) are characterized by 
direct coupling of the initial photons to the two different quark flavours. The complete set (II) 
is obtained from the representative diagrams shown in Fig. 1 by interchanging independently 
the photons and quark flavours. Of course, also set (II) is by itself gauge invariant. 




Fig. 1. Different topologies of the lowest- 
order Feynman diagrams contributing to 77 
Bcbc. 



More physically, the three subsets of diagrams discussed above can be interpreted as describing 
different production mechanisms. The subset (1;,) with primary 6-quarks may be associated 
with bb production and 6-quark fragmentation b B^c in accordance with the flrst diagram 
of Fig. 1 (I). Similarly, one may be tempted to attribute the subset (Ic) to cc production 
followed by c-quark fragmentation c — > Bcb. However, while the fragmentation description 
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is applicable to the process (I;,), it does not hold for (Ic) as will be shown later. In contrast 
to the class (I) diagrams, the class (II) diagrams are characterized by direct production of 
a bb and cc pair from the initial photons and recombination of the b- and c-quark into a Be 
meson, be —>■ Be. Obviously, the production amplitudes for the three different mechanisms are 
proportional to different combinations of the electric quark charges, namely Ql = 1/9, Ql = 
4/9 and QbQc = —2/9, respectively. As one of the main results, we will show later that the 
recombination mechanism (II) dominates, while bb radiation from a primary c-quark (Ic) gives 
a small contribution and fragmentation of a primary 6-quark (I^) is completely negligible. 

We begin with a brief description of technical details of the calculation of amplitudes and cross 
sections. For definiteness and because of their particular importance, we concentrate on the 
pseudoscalar and vector ground states. Be and B*, respectively. Because of the nonrelativistic 
nature of these systems, the relative momentum p of the b and c constituents and their binding 
energy are expected to be small in comparison to the quark masses. It is then reasonable to 
expand the amplitudes in p, and to keep only the lowest nonvanishing term 0. In this 
approximation, the b- and c-quark are on mass shell and move together with equal velocity. 
The following relations hold: 

TTl TTlh 

M{Be)=M{B:)=M = me + mb, Pc=j^P, Pl= J^P- (1) 

Moreover, the amplitudes for the production of S-waves reduce to hard scattering amplitudes 
multiplied by the S'-wave function \E'(0) at the origin. For P- waves one obtains a similar 
product involving the derivative of the P-wave function. Using the definitions 

< 0|67/.75c|Pe(P) > = ^fBP^. (2) 

and 

<Q\hAB:{P)>= MfB^^e, (3) 

one can substitute the nonrelativistic wave function parameter \E'(0) by the corresponding 
decay constants 

fB. = fBi = ^\^m- (4) 

As a simple rule, in the Feynman diagrams of Fig. 1 the transition from open quark production, 
77 — >■ bcbc, to B^*^ bound state production is achieved by making the following replacement 
1: 

/ ( ) 

v{p-,)u{pe) = ^{r-M)Uss,. (5) 

Here, v{pi) and u{pc) denote the Dirac spinors of the quarks forming the P^*-* bound states, 
and 1155^ is the appropriate spin projector, IIoo = 75 for pseudoscalars and IIi^^ = / for 
vectors. In Fig. 1, this formal substitution is indicated by the dark dot. Obviously, because 
of the neglect of the relative momentum of the constituent quarks, the Feynman amplitudes 
involve no loop in momentum space, but only traces of Dirac and colour matrices. Note that 
the colour structure is not accounted for in eq. (^. All diagrams of Fig. 1 have the same 
colour factor | which has been included in the overall normalization of the cross sections. 

The squared matrix element is obtained by two independent methods. One calculation is based 
on the traditional trace technique. The necessary traces are calculated using the symbolic 
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manipulation program FORM |Q. The resulting expressions are automatically implemented 
into a FORTRAN program. The second calculation is based on a generalization of the method 
described in detail in ref. |^ for bosonic final states in e~^e~ scattering. Here, the polarized 
matrix elements are calculated in a fixed Lorentz frame using the Weyl representation of the 
Dirac matrices and spinors together with real polarization vectors. The matrix element is 
reduced analytically to a linear combination of products of the 2x2 Pauli matrices. These 
products are computed numerically. The use of real photon polarization vectors allows to 
reduce the number of elements of the matrices which have to be calculated by a factor of two. 
The polarized matrix elements are obtained by sandwiching these matrices between 2x2 Pauli 
spinors. As the last step, one sums over polarizations. The matrix elements squared obtained 
by the two methods agree to ten digits. 

The phase space integration is also performed independently using two different routines. Be- 
fore employing the Monte Carlo integration routine VEGAS the strongest peaks of the matrix 
elements were eliminated by introducing new integration variables. After the integration all 
differential distributions and cross sections agree within the Monte Carlo errors. 



3 Results 

For numerical purposes and all subsequent figures, we choose the following values of the 
parameters: 

mb = 4.8GeV, = 1.5 GeV, a, = 0.2, a = 1/129, /^^ = /s* = 0.4 GeV [7]. (6) 

In Fig. 2 we compare the cross section a {lb + Ic + H) for 77 — > Bcbc derived from the complete 
set of O(a^af) diagrams indicated in Fig. 1 with the cross sections cr(Ife) and (j{lc) resulting 
from the gauge invariant subset (I) with primary b- and c-quarks, respectively. Also shown 
are the results obtained by simply multiplying the cross sections a^i and for open flavour 
production with the respective fragmentation probabilities P(b ^ Bcc) and P(c — BJ)) P . 
These probabilities have been calculated separately in the same nonrelativistic approximation 
as described in section 2. 

As can be seen, the total production cross section is dominated by the recombination mecha- 
nism except very near to the production threshold, where also fragmentation, or more precisely 
radiation, processes give significant contributions. Moreover, it is more efficient to produce 
a pair of c-quarks and to radiate from them a 66-pair in order to form a Be meson than 
to produce a pair of 6-quarks and to radiate a cc-pair. The reason is the electric charge of 
the primary quarks which yields a factor 16 in favour of c-quarks. Finally, the factorized 
description in terms of open flavour production, 77 — >■ hb and 77 cc, followed by the frag- 
mentation processes, h — > BcC and c Bcb, respectively, provides a good approximation to 
the class (I) processes only when the primary quark is a 6-quark, but not when it is a c-quark. 
These findings are in sharp contrast to B^ production in e"^e~-annihilation There, 
the recombination mechanism is absent and the radiation of a hh pair from primary c-quarks 
is negligible. Furthermore, the fragmentation picture works perfectly no matter whether the 
fragmenting quark is the heavier or the lighter one. 
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Fig. 2. Integrated cross sections for 77 
Bfbc versus the cm. energy. The different cal- 
culations are explained in the text. The produc- 
tion mechanisms are as classified in Fig. 1. 



We have investigated the break-down of the fragmentation description for the class (I) pro- 
cesses in 77-scattering in some more detail. In energy range considered, this approximation 
is indeed only applicable when the fragmenting quark is appreciably heavier than the quark 
flavour produced via gluon radiation. Such a favourable case is illustrated in Fig. 2 by (T{lb) 
and af,bP(b Bcc) which agree nicely. However, already for equally heavy quarks, a case 
realized in the analogous process 77 J/ipcc, we find considerable disagreement. To be 
more definite, (1(77 ^ .J /^cc) / acc2P{c J/ipc) ^ 2.1 (1.4) at = 20 (500) GeV. When 
the fragmenting quark is the lighter one, the fragmentation description fails completely as 
demonstrated by cr(Ic) and accP{c BJS) in Fig. 2. It is the second diagram of Fig. 1 (I) 
which seems to be responsible for the break-down of factorization. This diagram is absent in 
e^e~ annihilation, but it is necessary in 77-scattering for gauge invariance. 

Further features of the different production mechanisms and theoretical descriptions can be 
learned from differential distributions of the Be mesons. Fig. 3 shows the distributions in the 
scaled energy variable z = 2Eb^/ y/s obtained from the complete set of 0{a'^al) diagrams and 
from the gauge-invariant subsets (If,) and (Ic) alone. Most notably, the dominant recombina- 
tion processes populate the soft end of the spectrum, while the class (I) processes give rise 
to a characteristic high energy tail. The hard component (lb) is as expected from 6-quark 
fragmentation. In fact, the energy distribution predicted from the diagrams of class (I;,) is 
practically indistinguishable from the corresponding fragmentation function Di^^^ {z) ||, . 
However, for the component (Ic) the analogous distributions do not coincide. The actual en- 
ergy spectrum derived from the diagrams (Ic) is considerably harder than the fragmentation 
function Dc^bX^) B I13- '^^^ above assertions are illustrated more clearly in Fig. 4 where we 



plot the individual energy spectra generated by the processes (If,) , (Ic) and (II) in comparison 
to the fragmentation functions Di^gX^) and Dc^bX^)- Here, all distributions are normalized 
to unity. 
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Fig. 3. Energy distributions of Be mesons at 
^/s = 100 GeV with all processes of Fig. 1 in- 
cluded and from the radiative mechanisms (If^) 
and (Ic) alone. 



Fig. 4. Normalized energy spectra of Be 
mesons at y/s = 100 GeV generated by the 
mechanisms (II), (Ih) and (Ic) of Fig. 1 in com- 
parison to the heavy quark fragmentation func- 
tions. 



Interesting observations can also be made from the pt distributions depicted in Fig. 5. The 
factorized description in terms of hh production and 6-fragmentation and the corresponding 
diagrams (It) of Fig. 1 yield distributions which are very similar in shape and normalization 
except in the low-p^ region. As far as the shape is concerned this is even true for the analogous 
processes with primary c-quarks. However, the c-quark fragmentation picture fails to repro- 
duce the correct magnitude of the cross section for the process (Ic) , a problem already noticed 
in Fig. 2. We also see that the radiation of a cc-pair from a 6-quark leads to a harder pt 
distribution for the bound states than the radiation of a 66-pair from c-quarks. Finally, the 
most important observation is that the recombination mechanism dominates production 
not only at low pr, as one could have expected, but also in the high-p^ region up to the kine- 
matic limit. In other words, the familiar description of high-p^ hadron production in terms of 
the production and fragmentation of quarks with large transverse momenta is inadequate for 
single B(. production in 77-scattering and similar processes. 

For completeness, we mention without showing figures that the angular distributions resulting 
from the different mechanisms considered above are all quite similar in shape. The distribu- 
tions are sharply peaked in forward and backward direction, with the distribution predicted 
by the complete set of diagrams of Fig. 1 being slightly flatter than the distributions derived 
from subsets (It) and (Ic), and from h- and c-quark fragmentation. 

So far we have concentrated on the problem of identifying the dominant production mech- 
anisms and examining different theoretical descriptions. For this purpose, it is sufficient to 
consider photon-photon collisions at fixed energies. Now we want to make predictions on the 
production rates of Be mesons in photon-photon collisions at e"*"e~-machines in the LEP energy 
range and beyond. To this end we fold the total cross sections for 77 — > B^^hc with realistic 
photon spectra that is the photon spectrum obtained by Compton back-scattering of high 
intensity laser light on beams or the Weizsacker- Williams bremsstrahlung spectrum 
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[0. Our results are shown in Fig. 6, where we have plotted the convoluted cross sections for 
production of the pseudoscalar and vector ground states versus the e^e~ centre-of-mass en- 
ergy. For reference we also show the unfolded cross sections for the subprocesses 77 B^*^bc 
as functions of the 77 centre-of-mass energy. As expected from the shape of the Compton 
spectrum which has a long tail towards soft photons, and the shape of the 77-cross section 
which peaks just above threshold, the convolution increases the cross sections substantially for 
energies above 100 GeV. At a 500 GeV linear collider and for an integrated luminosity of 10 
fb~^, one can produce about 100 Be and 400 B* mesons. In contrast, the yield of Be mesons 
from bremsstrahlung photons is very small at LEP energies, but increases logarithmically with 
energy. While at LEP 2 the production rate is still too low to lead to observable signals, in 
the TeV energy range bremsstrahlung photons become competitive with back-scattered laser 
photons in producing B^ mesons. 




Fig. 5. The Be transverse momentum dis- 
tributions at ^/s = 100 GeV for the differ- 
ent mechanisms and descriptions considered in 
Fig. 2. 



Fig. 6. Total production cross sections for B* 
(upper curves) and Be (lower curves) versus 
the 77 cm. energy (full), and convoluted with 
the Compton spectra of back- scattered laser 
light (long- dashed) and with bremsstrahlung 
spectra (short- dashed) versus the e+e~ cm. 
energy. 



4 Summary 

We have studied the production of Be mesons in photon-photon scattering. Our main results 
can be summarized as follows. The dominant production mechanism proceeds via direct 
coupling of a bb and cc pair to the initial photons and recombination of the b- and c-quark 
into the Be bound state. The radiative process where a cc-pair is produced via gluon emission 
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from a primary 6-quark plays a negligible role. It is well described by 6-quark production and 
fragmentation, 77 — > 66, 6 — BcC. The analogous process where a 66-pair is emitted from 
a primary c-quark is enhanced by the charge ratio Q\/Ql, but it still plays a marginal role. 
Here, the fragmentation description 77 — > cc, c ^ Bch does not hold. 

Given the dominance of recombination even at large transverse momenta up to the kine- 
matic limit we conclude that the usual hard scattering formalism in terms of high-p^- quark 
production and subsequent fragmentation cannot be applied to 77 B^X. This also puts 
serious doubts on the calculations of hadronic B^ |]13| and J/ ip [1^ production based on bb 



and cc production via gluon-gluon fusion followed by heavy quark fragmentation. Although 
the replacement of the electric charges by colour charges when going from 77- to (7 (7-pro cesses 
changes the relative weights of the production mechanisms considered in this paper, we sus- 
pect the fragmentation description to underestimate the yields of heavy bound states in the 
experimentally accessible, high-p^ regime. 

The Bj^*^ production rates predicted for 77 collisions and realistic photon energy distributions 
show that these processes are below the observability limit in e"'"e~-collisions at LEP energies. 
However, at linear colliders in the TeV energy range, two-photon production of B^ mesons 
may come into experimental reach. These general prospects are little affected by the sizeable 
uncertainties in the values of parameters quoted in eq. (^, and the unknown higher order 
perturbative, nonperturbative and relativistic corrections. 
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